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ABSTRACT 



We present a detailed model for the ionized absorbing gas evident in the 900 
^ I ksec Chandra HETGS spectrum of NGC 3783. The analysis was carried out 

Q ' with PHASE, a new tool designed to model X-ray and UV absorption features 

^ ■ in ionized plasmas. The 0.5-10 keV intrinsic continuum of the source is well 

^ . represented by a single power law (F = 1.53) and a soft black-body component 

O ! (kT ~ 0.1 keV). The spectrum contains over 100 features, which are well fit by 

g I PHASE with just six free parameters. The model consists of a simple two phase 

^ ' absorber with difference of ~ 35 in the ionization parameter and difference of 

Q^' 4 in the column density of the phases. The two absorption components turned 

Q . out to be in pressure equilibrium, and are consistent with a single outflow 750 

^ ! km s~^) and a single turbulent velocity (300 km s~^), and with solar elemental 

^ I abundances. The main features of the low ionization phase are an Fe M-shell 

unresolved transition array (UTA) and the O Vll lines. The O Vll features, 
usually identified with the O vill and a warm absorber, are instead produced 
in a cooler medium also producing O vi lines. The UTA sets tight constraints 
on the ionization degree of the absorbers, making the model more reliable. The 
high ionization phase is required by the O vill and the Fe L-shell lines, and there 
is evidence for an even more ionized component in the spectrum. A continuous 
range of ionization parameters is disfavored by the fits, particularly to the UTA. 
Our model indicates a severe blending of the absorption and emission lines, as well 
as strong saturation of the most intense O absorption lines. This is in agreement 
with the O Vll {t\ = 0.33) and O vill (ta = 0.13) absorption edges required 
to fit the spectrum. The low ionization phase can be decomposed into three 
subcomponents, based on the outflow velocity, FWHM, and H column densities 
found for three out of the four UV absorbers detected in NGC 3783. However, the 
ionization parameters are systematically smaller in our model than derived from 
UV data, indicating a lower degree of ionization. Finally, our model predicts a Ca 
XVI line for the feature observed at around 21.6 A (a feature formerly identified 
as O vil) , constraining the contribution from a zero redshift absorber. 



Subject headings: galaxies: absorption linjs - galaxies: Seyferts - galaxies: active - 
galaxies: X-ray 



1. Introduction 

Quasars and Active Galactic Nuclei 
(AGNs) display a rich and confusing array 
of atomic emission and absorption features 
throughout the infrared, optical, ultravio- 
let and soft X-ray bands. Does this pro- 
fusion of features reflect an inherent com- 
plexity or even randomness within AGNs? 
Or is there an underlying order that is 
masked by the overabundance of detail? 
The case has been made on both sides (for 
complexity: Baldwin et al. 1996; Korista 
et al. 1996; Krohk & Kriss 2001; for under- 
lying order: Mathur, Elvis & Wilkes 1995; 
Elvis 2000; Arav et al. 2003). 

One case has been intensively pursued 
lately, due to the improvement in the 
quality of UV and X-ray spectra from 
STIS, FUSE, and the grating spectrome- 
ters on Chandra (HETGS, LETGS) and 
XMM-Newton (RGS): the highly ion- 
ized (or 'warm') absorbers that appear 
in about half of all UV and X-ray spec- 
tra of AGNs. Both UV and X-ray ab- 
sorption lines show moderate velocity 
(~1000 km s~^) blueshifts, implying out- 
flows. The mass loss rate from these out- 
flows can be a substantial fraction of the 
accretion rate needed to power the AGN 
continuum, and so these 'quasar winds' 
are likely to be dynamically important to 
AGNs (Mathur, Elvis & Wilkes 1995). Un- 
derstanding the physical state or states of 
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this gas, and hence the location and dy- 
namics of AGN winds would be a valuable 
advance in the understanding of quasars. 

The UV and X-ray absorbers must be 
closely related: exactly the same 50% of 
AGNs show both (Mathur, Elvis, & Wilkes 
1995; Mathur, Wilkes, & Elvis 1998; Cren- 
shaw et al. 1999; Monier et al. 2001), 
and each can be used to predict the pres- 
ence of the other. However, the absorbers 
show complex, multi-component and time- 
variable structures in UV spectra, and it 
is not straightforward to connect these to 
the X-ray lines, given the order of magni- 
tude lower spectral resolution in the X-ray 
band. Moreover column densities derived 
from the UV lines are strongly affected by 
scattered light into our beam or by partial 
covering of the continuum source (Arav et 
al. 2003). While the UV absorbers are 
dominated by O VI hues, with C IV and N 
V also appearing, the X-ray absorbers are 
dominated by O VII and O VIII, with many 
other lines of similar ionization now be- 
ing seen. Despite being dubbed 'warm ab- 
sorbers' the absorbers respond to changes 
in the ionizing continuum and so are at 
least in part photoionized (Nicastro, Fiore, 
& Matt 1999). Hence we shall use the less 
loaded term 'ionized absorbers'. 

Here we take the best existing high res- 
olution X-ray spectrum of an ionized ab- 
sorber in an AGN, a 900 kscc observation 
of NGC 3783 with the Chandra HETGS, 
and analyze the absorption features self- 
consistently with a new code ('PHASE', 
Krongold et al. 2003, in preparation). We 
flnd a strikingly simple solution: well over 
100 absorption features can be fitted by 
a two component solution having only six 
free parameters. Moreover, allowing for 
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the lower X-ray resolution, the lower ion- 
ization component is consistent with two 
UV absorbers. There is also evidence in 
the X-rays for a third, high velocity and 
low column density, UV component in four 
lines. We discuss in detail why we are 
able to obtain this simple solution, in con- 
trast to earlier studies of the same object 
and spectrum that led to more complex 
answers. The two component solution is 
strongly suggestive of a two-phase medium 
in pressure balance, with further likely im- 
plications. 

1.1. High Resolution X-ray Spectra 
of NGC 3783 

NGC 3783 is an extensively studied 
bright Seyfert galaxy, at a redshift 0.00976 
±0.00009 (2926 ±28 km s-^ de Vau- 
couleurs et al. 1991). NGC 3783 has been 
the subject of extensive monitoring by the 
High Energy Transmission Grating Spec- 
trometer (HETGS, Canizares et al. 2000) 
on board the Chandra X-ray Observatory, 
with a total exposure of ~ 900 ksec (re- 
ported in Kaspi et al. 2002). The result- 
ing spectrum covers the 0.5-10 keV energy 
range, and is the best obtained so far for 
an ionized absorber in a Seyfert galaxy, 
with more than 2000 counts per resolution 
element at 7 A. This spectrum showed the 
presence of more than one hundred ab- 
sorption lines from a wide range of ionized 
species. A few weak emission lines, mainly 
from O and N, were also detected. The 
spectrum further revealed a wide absorp- 
tion feature in the 16 — 17 A range. Kaspi 
et al. (2002) identified this feature as an 
unresolved transition array (UTA), aris- 
ing from numerous Fe M-shell inner-shell 
2p — 3d transitions. The XMM-Newton 



RGS has also observed NGC 3783 (Blustin 

et al. 2002). This observation extends to 
longer wavelengths than the 900 ksec expo- 
sure from the HETGS. A complete review 
of the analyses carried out on these data 
will be presented in §5.5. To this date, self 
consistent modelling for the full 900 ksec 
Chandra high resolution spectrum has not 
been published. 

NGC 3783 has also been studied in 
the UV band with FUSE and HST-STIS. 
These observations have revealed intrin- 
sic absorption by Lya, O VI, and N V 
in this region of the spectrum (Kraemer, 
Crenshaw & Gabel, 2001, and references 
therein). The high spectral resolution and 
signal to noise ratio of both data sets al- 
lowed the identification of four absorption 
systems moving toward us, with outfiow 
velocities of 1320, 1027, 724, and 548 km 
s-^ (Gabel et al. 2003). 

In this paper, we present the first self- 
consistent model for the 900 ksec X-ray ab- 
sorption spectrum of NGC 3783. In the 
following section (§2) we describe the data 
reduction. The modelling was carried out 
with "PHASE", a newly developed code 
based on CLOUDY (v. 90.04, Ferland 
1997) and the Astrophysical Plasma Emis- 
sion Database (APED v. 1.3.0, Smith et al. 
2001), and designed to reproduce absorp- 
tion features produced by ionized plasmas. 
This code is briefiy explained in section §3. 
In §4 we present the results obtained for 
the spectral fitting and in §5 and §6 we 
discuss the rehabifity of the model and its 
implications. 

2. Data Reduction 

NGC 3783 was observed six times us- 
ing the HETGS on board the Chandra X- 
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Fig. 1. — Light curve of Chandra-HETG 
observation of NGC3783 in the 2.5-4 keV 
(3-5 A) range. 

ray Observatory, with the Advanced CCD 
Imaging Spectrometer (ACIS). The first 
observation was carried out on January 
2000 and had a duration of 56 ksec (re- 
ported in Kaspi et ah 2000). Between 
February and June of 2001, the following 
five observations were taken, each one with 
an exposure time of 170 ksec (reported in 
Kaspi et al. 2002). Figure 1 presents 
the light curve in the 2.5-4 keV (~ 3.0- 
5.0 A) range. We chose this range because 
it avoids the Fe K 6.4 keV complex and 
most absorption from bound-free transi- 
tions. Therefore, it can indicate changes 
from the central source. In agreement with 
Kaspi et al. (2002), we found variability at 
a level of about 50% between the fourth 
and fifth observations, and smaller in the 
rest of the cases. In this paper we present a 
model for the averaged spectrum of the six 
observations. We do not expect the level of 



variability to affect our main results. How- 
ever, a detailed analysis on each observa- 
tion, and the general effects of variability 
will be presented in the future. 

We retrieved from the public Chandra 
data archive^ the primary and secondary 
data products for these six observations 
and reprocessed their event files with the 
Chandra Interactive Analysis of Observa- 
tions (CIAO^) software (Version 2.3). Fol- 
lowing the on-line data analysis "threads" , 
provided by the Chandra X-ray Center 
(CXC^), we extracted from the CALDB^ 
(Gracsslc et al. 1996) the source spectra 
and the corresponding Ist-order Ancillary 
Response Files (ARFs) and Redistribution 
Matrices (RMFs). 

The HETGS (Canizares et al. 2000) 
produces spectra from two grating assem- 
blies, the high-energy (HEG) and medium 
energy (MEG) gratings. All the spectra 
described below and throughout the paper 
correspond to the MEG. We combined the 
-1st and -|-lst orders for each exposure. To 
further increase the signal to noise ratio of 
the data, we coadded the spectra and av- 
eraged the ARFs, weighting them accord- 
ing to their respective exposure times. The 
resulting spectrum has a net exposure of 
888.7 ksec. The signal to noise ratio per 
resolution element (^^0.02 A for the MEG) 
varies from 46 to ~ 5 in the 6-25A 
0.5-2.0 keV) range. The fluxed spectrum is 
presented in Figure 2a, and the empirical 
spectrum is presented in Figure 3. 



^http 
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/ / asc.harvard.edu/ cda 

//asc. harvard.edu/ciao 

//asc. harvard.edu/ciao/documents_t hreads.html 
/ / cxc.harvard.edu/ caldb 
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Fig. 2. — Fluxed MEG first-order spec- 
trum of NGC 3783. (a) Continuum fit- 
ted to the spectrum. The three sohd hues 
stand for (from upper to lower): intrinsic 
continuum of the source, continuum atten- 
uated by Galactic absorption, and contin- 
uum further attenuated by the two phase 
absorber edges. The dotted line represents 
the predicted power law without the con- 
tribution from the black body component, 
(b) Two phase absorption model plotted 
for comparison.^ 

^Figures available in color in the electronic format. 



3. Modelling: PHASE 

With the aim of modelling in a self- 
consistent way all the absorption features 
observed in the spectra of UV and X-ray 
sources, we have developed PHASE (PHo- 
toionized Absorption Spectral Engine), a 
code that calculates absorption due to 
an ionized plasma based on ATOMDB^ 



(v.1.3.0. Smith et al. 2001) and CLOUDY 
(v. 90.04, Ferland 1997). The calcula- 
tions assume a simple geometry that con- 
sists of a central source emitting an ion- 
izing continuum with clouds of gas in- 
tercepting our line of sight. The ioniza- 
tion balance is calculated using CLOUDY 
with an SED chosen to be consistent with 
the intrinsic Spectral Energy Distribution 
(SED) of the source. The parameters of 
the code are (1) the ionization parame- 
ter (defined over the entire Lyman con- 
tinuum, as [/ = Q{H)/47Tr'^n{H)c with 
Q{H) the number of H ionizing photons 
per second, r the distance to the source, 
n{H) the H density, and c the speed of 
light.), (2) the equivalent hydrogen col- 
umn density, (3) the outflow velocity, and 

(4) the internal micro-turbulent velocity 
of each one of the absorbers, as well as 

(5) the intrinsic SED of the source. A 
sixth parameter, the electron temperature, 
can be included, forcing the temperature 
to be different from photoionization equi- 
librium. PHASE can thus handle hybrid 
states, e.g. gas photoionized by a strong 
radiative source but where another exter- 
nal mechanism (for instance, shock heat- 
ing) keep the electron temperature at a 
higher value than the one expected only 
through photoionization equilibrium. The 
micro-turbulent velocity and the thermal 
velocity are added in quadrature to com- 
pute the total Doppler broadening param- 
eter of the gas: V^op = V^^^^^ + l/^^^^. 

The code proceeds as follows: (1) The 
absorption produced by each transition of 
the different ion species is calculated, using 



^http: / / cxc.harvard.edu/ atomdb 
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the formula 

Jo rnc 

where ttj, is the absorption coefficient at 
the frequency u, L is the hnear size of 
the cloud, Nion is the column density of 
the given ion (calculated through the ion- 
ization balance), m is the electron mass, 
e is the electron charge, c is the speed 
of light, flu is the oscillator strength of 
the electron transition (from the lower to 
the upper level), and is the line proffie 
(which is assumed to be a Voigt profile, 
Rybicki & Lightman, 1979, and calculated 
by convolving the natural broadening with 
the thermal and microturbulent Doppler 
broadening) . 

(2) Then the model calculates the con- 
tinuum absorption produced by bound-free 
transitions. 

(3) The contributions from different ab- 
sorption lines and edges are added together 
to produce an absorption spectrum. 

(4) This spectrum is finally shifted to 
the proper outflow velocity. 

PHASE includes ~ 4000 lines arising 
from the ground level. The line list was 
extracted from the Astrophysical Plasma 
Emission Database (APED v. 1.3.0, Smith 
et al. 2001). APED includes informa- 
tion up to the 5th principal quantum 
number for transitions arising from the 
14 more abundant elements in the uni- 
verse, namely H, He, C, N, O, Ne, Mg, Al, 
Si, S, Ar, Ca, Fe, and Ni (for a detailed 
list of the ions of each element included 
in APED visit the ATOMDB web page, 
http://cxc.harvard.edu/atomdb). To fur- 
ther increase the number of lines in our 
model, we added nearly 200 features in- 



volving principal quantum numbers n=5 
to n=10 (levels not present in APED) 
from the "Atomic Data for Permitted Res- 
onance Lines" (Verner, Verner, & Fer- 
land 1996) . This compilation includes 
lines connecting the ground level to ex- 
cited levels via optically allowed transi- 
tions, i.e. transitions between levels of 
opposite parity with |AL| < 1, [A^*! = 0, 
and |AJ| < 1, where L is the orbital, S 
the spin, and J the total angular momen- 
tum quantum number. In order to include 
in our calculations the Fe M-shell inner 
transitions (which are responsible for the 
UTA), we used the abbreviated data set 
provided by Behar, Sako, Sz Kahn (2001), 
the best source of information for these 
transitions available in the literature. This 
abbreviated version of the data assumes 
that mean wavelengths and effective oscil- 
lator strengths can emulate the complex 
absorption produced by several lines of 
each ion. According to the authors, the 
approximation can successfully reproduce 
the general shape and EW of the feature. 
The only exception is for Fe x-xill ions; 
in this case, if the ion column density lies 
in the range 10^'' — 10^^ cm~^ and the tur- 
bulent velocity is < 100 km s~^, the ap- 
proximation can overestimate the actual 
value of the EW by as much as 40%. We 
also included the most intense (oscillator 
strengths > 0.1 mA) inner shell {ls-2p) 
transitions from Ne, Mg, Al, Si, S, Ar, Ca, 
and Fe. These data were extracted from 
the list of Behar & Netzer (2002), and will 
be included in APED (v.2.0.0). A fuh ex- 
planation of the present code will be pub- 
lished in a forthcoming paper (Krongold et 
al., in preparation). 

We performed spectral fitting of the 
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NGC 3783 data using PHASE integrated 

into the Sherpa (Freeman et al. 2001) 
package in CIAO (Fruscione2002) as a ta- 
ble model. 

4. Spectral Fitting 

The key point of our modelling was to fit 
self-consistently the absorption hues and 
edges from all the ions of all the species 
at the same time, without any assumption 
about the absorption of individual ions. 
Therefore, we did not introduce any prior 
constraint on the column density or popu- 
lation fraction of any ion. As will be dis- 
cussed in §5, many lines of different ions 
are blended, making the measurement of 
individual ion column densities unreliable 
in most cases. The determination of in- 
dividual ion column densities through con- 
tinuum bound-free absorption is unreliable 
also due to the contribution of so many dif- 
ferent ions. Furthermore, even if the col- 
umn densities produce edges, these might 
be hidden by absorption lines. Such a case 
is clearly exemplified in Figure 4 for the 
spectrum of NGC 3783: the O VII edge 
and the UTA at 16 — 17 A are completely 
blended, making very difficult to isolate 
the contribution of the former to the ab- 
sorption. Our model is self-consistent in 
the sense that given the intrinsic SED of 
the source, and the column density and 
ionization parameter of the absorbing me- 
dia, the ionization balance is fixed. By 
also including the broadening mechanisms 
of the line as well as the rules governing 
the ion transitions, in a self-consistent way, 
a final global solution consistent with all 
these processes, can be obtained. This pro- 
cedure avoids the inconsistencies that may 
arise in the ionization structure of the gas. 



as well as the effects of neglecting satu- 
ration, that are present in analyses where 
only an ion by ion (or line by line) approach 
is used to fit the spectral features. 

The intrinsic continuum and the gas ab- 
sorption were fitted simultaneously. There- 
fore, we iterated many times to get the 
final X-ray SED. We further attenuated 
the continuum by an equivalent hydro- 
gen column density of 1.013 x 10^^ cm~^ 
(Murphy et al. 1996) to account for the 
Galactic absorption by cold gas in the di- 
rection to the source. In the calculations 
reported here, we have explored only pho- 
toionization equilibrium models, and we 
have assumed solar elemental abundances 
(Grevesse, Noels, & Sauval 1993). 

4.1. Continuum 

To model the intrinsic continuum of the 
source, we first fitted a simple power law 
with varying photon index and amplitude. 
However, a single power law could not 
fit the data over the entire range. We 
found an excess of fiux in the spectrum 
at energies < 0.6 keV. This excess could 
be accounted for with the inclusion of a 
black body component superimposed on 
the power law. Such a thermal component 
has been used in the past to describe the 
NGC 3783 continuum (e.g. De Rosa et al. 
2002, and references therein). 

To further check our results, we made 
an independent estimate of the continuum 
level using the wavelength bands where 
no line features are present, the "line- free 
spectral bands" (horizontal green ranges 
at the bottom of Fig. 3 panels). These 
bands were determined excluding the re- 
gions with lines in our model, plus the ob- 
vious features that our model missed 
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Fig. 3. — Two phase absorber model plotted against the first-order MEG spectrum of NGC 
3783. Absorption lines predicted are marked in the top (red). Single labels stand for emission 
lines (blue). The line free zones are indicated at the bottom of each panel (green). The 
continuum level (including edge continuum absorption) is overplotted for comparison (dotted 
green line). The spectrum is presented in the rest frame system of the absorbing gas.^ 
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(see §5). The two continuum estimates 
were consistent with each other. Table 1 
shows our resuhs. Figure 2a presents the 
intrinsic continuum of the source, the con- 
tinuum attenuated by Galactic absorption, 
and the observed continuum (i.e. the con- 
tinuum further attenuated by the bound- 
free transitions from our model). As shown 
in this figure (and also in Fig. 3) there is 
an excess of flux in the observed continuum 
at around 5.5-6.5 A. We attribute this ex- 
cess to a calibration effect due to the rapid 
changes of the effective area between 1.7 
and 2.3 keV, resulting from the Iridium M- 
edge. As can be inferred from the figure, a 
lower photon index would account for the 
difference, but would overpredict the con- 
tinuum level at shorter wavelengths. 

For simplicity, to calculate the ioniza- 
tion balance, we approximated the full 
SED as a power law with the following 
photon indices: F = 2 below 13.6 eV, 
F = 2.4 in the range 13.6 eV < hu < 
0.6 keV, and F = 1.53 above 0.6 keV (i.e. 
we matched the UV SED derived from UV 
studies, see below). We further restricted 
our SED including high and low energy 
cutoffs (F = 5 above 100 keV and F = —3.5 
below 0.1 eV). The SED is plotted in Fig- 
ure 5. With the chosen continuum, be- 
tween 0.6 keV and 0.1 eV, we assumed the 
same SED used in Kraemer, Crenshaw, & 
Gabel (2001) for their modelling of the UV 
spectrum of NGC 3783. Although this ap- 
proximation appears reasonable, there are 
important uncertainties associated with it, 
due to the unobservable region in the UV- 
X-ray range associated with Galactic ab- 
sorption and the lack of simultaneous ob- 
servations in both bands. The shape of 
the SED in the UV and far UV does not 



affect the absorption in the X-ray range, as 
discussed by Kaspi et al. (2001) and Steen- 
brugge et al. (2003). However, the exact 
shape of the SED does affect the values 
predicted for the ionization parameter, as 
will be discussed in §5.4. 

4.2. Emission Features 

Kaspi et al. (2002) identified several 
emission lines apparent in the 900 ksec 
spectrum. In our model, these emission 
lines partly blend with the absorption 
lines, such that the depths and EWs of 
both appear reduced when the absorption 
and emission components are added to- 
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Fig. 4.— PHASE model for NGC 3783 
warm absorber. The top dotted line (red) 
shows the O K-edges due to bound free ab- 
sorption. The lower solid curve includes 
additional contribution from bound-bound 
transitions. As can be seen, without the 
inclusion of the resonant absorption, the 
optical depth of the edges can be severely 
overpredicted.^ 
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Table 1 
Continuum Parameters 



Power Law 

Photon Index (F) Normalization'^ NhguI (cm~^) 
1.53±0.02 0.011±0.002 1.013 x lO^^ 

Thermal component 

kT (keV) Normalization^ 
O.lOiO.03 0.00020±0.00007 



^In photons keV ^ cm ^ s ^ at 1 keV. 

^In Lsg/D^Q, where L39 is the source luminosity in 
units of 10^^ erg and Diq is the distance to the source 
in units of 10 kpc. 



gether (see Fig. 6). In its present stage, 
our code cannot yet model emission in a 
self-consistent way. To do so requires mod- 
elling the geometry and kinematics of the 
absorbing gas. Therefore, to account for 
the filling effect of the absorption lines, 
we fitted gaussian profiles to the 14 most 
prominent emission features in the spec- 
trum, and added them to the model. To 
make the fitting as constrained as possible 
we used the same outflow velocities and 
the same FWHMs to model all the emis- 
sion lines produced by the same ion. The 
emission lines fitted and their properties 
are listed in Table 2. We excluded the Fe 
K 6.4 keV complex. 

Two features in the spectrum have been 
tentatively identified as O Vll and N vi 
radiative recombination continua (RRC), 
because of their wavelengths. By contrast 
to hot coUisionally ionized plasmas where 



the RRCs are broad extended features, in 
photoionized plasmas the RRCs are narrow 
"linehke" features. The "width" of these 
features can be used, in principle, to esti- 
mate the electron temperature (Liedahl & 
Paerels 1996). However, due to the heavy 
blending of these features with absorption 
components, a measurement of the tem- 
perature from the RRCs would be unre- 
liable. To account for this effect, a self- 
consistent model for emission processes is 
also required. 

Apart from these two features, all 
the emission lines have outflow velocities 
smaller than those of the absorption com- 
ponents, and are consistent with being at 
rest in the system frame (in accordance 
with Kaspi et al. 2002). This points to a 
symmetric distribution of the emitting gas. 
Using the O Vll triplet plasma diagnos- 
tics by Porquet & Dubau (2000), an upper 
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Fig. 5. — Spectral Energy Distribution 
used to approximate the NGC 3783 con- 
tinuum (solid line). The SED by Mathews 
& Ferland (1987) is also plotted for com- 
parison (dashed line).^ 

limit for the electron density in the O VII 
line formation region of Ug < 2.5 x 10^° 
cm~^ can be imposed, consistent with the 
one obtained by Morales et al. (2003, in 
preparation, Ug < 2 x 10^° cm~'^). 

4.3. Absorption Features 

4-3.1. Components outflowing at 750 km 

As was the case for previous studies, 
the spectrum of NGC 3783 could not be 
reproduced by a single absorbing compo- 
nent. Therefore, we included a second 
one with different ionization. Both com- 
ponents were consistent with an outflow 
velocity of ~750 km s~^. Since the ab- 
sorption lines in the spectrum are unre- 
solved by the HETGS and many of them 



Fig. 6. — Absorption/emission blending 
for the O Vlll(A18.969A) line. While the 
dotted line stands for the emission com- 
ponent, the dashed one represents absorp- 
tion. The final result (our model) is plotted 
with the solid line. As can be observed, the 
measurement of EWs from the data (e.g. 
Kaspi et al. 2002) underestimate both ab- 
sorption and emission.^ 

are blended, it was impossible to deter- 
mine the turbulent velocity of the systems. 
Therefore, in order to fit the data, we set 
this velocity equal to 300 km s~^, the same 
value used in Kaspi et al. (2001). The tem- 
perature in each gas component was calcu- 
lated from our models assuming photoion- 
ization equilibrium, and is ~ 9.5 x 10^ K 
for the hotter component and ~ 2.6 x 10^ K 
for the cooler one; therefore, the contribu- 
tion of the thermal velocity to the Doppler 
velocity is rather small (< 10 % for oxy- 
gen). Our results are listed in Table 3. The 
plot with the model is presented in Figure 
3 (and also in Fig. 2b, where we present 
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Table 2 
Emission Lines Parameters 



Ion 


Rest Frame 


FWHM 


Outflow Vel. 


EW» 




A(A) 


(km s"-*-) 


(km s~^) 


(mA) 



Neix 


13.699 


406±169 


267±110 


15.3 


Oviii 


16.006 


542±343 


102±148 


25.3 


RRC Ovii 


16.771 


715±277 


1370±173 


16.2 


Ovii 


17.200 


632±269 


35±232 


7.3 


Ovii 


17.396 


632±203 


35±272 


30.3 


Ovii 


17.768 


632±298 


35±150 


34.3 


Ovii 


18.627 


632±222 


35±148 


25.8 


Oviii 


18.969 


542±235 


102±79 


79.0 


Nvii 


20.910 


370d=248 


198±130 


23.7 


Ovii 


21.602 


632±417 


35±208 


115.1 


Ovii 


21.807 


632±267 


35±265 


60.0 


Ovii 


22.101 


632±258 


35±124 


144.9 


RRC Nvi 


22.458 


333±302 


1030±106 


49.8 


Nvii 


24.781 


370±281 


198±387 


43.2 



+14.6 

-9.8 

+18.4 

-16.9 

+17.6 

-14.9 

+8.3 

-8.8 

+12.6 

-15.2 

+23.4 

-14.6 

+10.2 

-17.9 

+41.8 

-34.7 

+33.3 

-27.1 

+47.4.1 

-37.5 

+42.5 

-31.6 

+59.1 

-68.8 

+19.4 

-25.4 

+37.8 

-41.2 



^Uncertain due to blending with absorption lines 
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Fig. 7. — Absorption lines produced by the 
high ionization phase (HIP) of our model, 
plotted against the first-order MEG spec- 
trum of NGC 3783. The continuum level 
contains edge absorption from both com- 
ponents. As can be observed, below 16 A 
the spectrum is dominated by Fe L-shell 
absorption lines. ^ 

the model for the full spectral range plot- 
ted over the fluxed spectrum). 

We find no evidence for cosmic dust or 
anomalous abundances in any of the ab- 
sorbers. It is noteworthy that the outflow 
of both components is consistent with a 
single velocity. While kinematically indis- 
tinguishable and with column densities a 
factor 4 apart, our two components have 
dramatically different ionization parame- 
ters: the high ionization phase (HIP) has 
a value ~ 35 times larger than the low ion- 
ization phase (LIP). The HIP, presented in 
Figure 7, gives rise to the absorption by N 
VII, O VIII, Ne ix-x, Mg x-xii, Al xii- 
XIII, Si xi-xiv, S xiii-xvi, Ar xii-xvii, Ca 



Fig. 8. — Absorption lines produced by the 
low ionization phase (LIP) of our model, 
plotted against the first-order MEG spec- 
trum of NGC 3783. The continuum level 
contains edge absorption from both com- 
ponents. The main signature of this phase 
is the UTA, as well as the O vil lines. ^ 

xiii-xviii, Fe xvii-xxii, and Ni xix-xxi. 
On the other hand, the LIP (presented in 
Fig. 8) produces absorption by C vi, N 
VI, O VII, Ne v-viii, Mg vii-ix. Si vii-ix, 
and the Fe M-shell inner transitions form- 
ing the UTA (Fe vil-xil). The most in- 
tense lines, as well as the unblended fea- 
tures, are listed in Table 4, where we show 
the contribution from the HIP and the LIP 
separately. All the EWs derived from our 
model for unblended lines agree with those 
measured from the data (see Fig. 9). 

4-3.2. Component outflowing at 1345 km 

In the 900 ksec spectrum Kaspi et al. 
(2002) found evidence for two X-ray ab- 



14 



Table 3 

Two Phase Absorber Parameters 



Parameter High-Ionization Low-Ionization 



Log U'^ 


0.76±0.1 


-0.78±0.13 


Log (cm-2)" 


22.20±.22 


21.61±0.14 


Vrurb (km s"^) 


300 


300 


Yout (km s-^r 


788±138 


750±138 


T (K)^ 


9.52 ± 0.44 X 10^ 


2.58 ±0.39 X 10^ 


[Log T (K)] 


5.98±0.02 


4.41±0.07 


Log T/U (ocP'^) 
Log 


5.22±0.12 


5.19±0.20 


1.02±0.12 


0.99±0.20 



'^Free parameters of the model. 

^^Derived from the column density and ionization pa- 
rameter, assuming photoionization equilibrium. 

'^The pressure PocUgT. Assuming that both phases lie at 
the same distance from the central source ocl/U, and 
PocT/U. 

^Ionization Parameter as in Krolik & Kallman 1984), 
i.e. the ratio between the radiation pressure and the gas 
pressure. 
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Fig. 9. — Measured EWs vs. model pre- 
dicted EWs for unblended lines in the NGC 
3783 spectrum. A line with slope equal to 
one is presented for comparison. 

sorbers in the hues of O Vll, but could 
not conclude the same in those of Ne X. 
We also see evidence for these O VII lines 
(lower panels of Fig. 10). As can be in- 
ferred from Table 4, while our model indi- 
cates that all the absorption produced by 
Ne X arises in the HIP, most of the O Vll 
absorption originates in the LIP. Further- 
more, no evidence of a second absorber in 
O VIII is found, and this ion too is mainly 
produced by the HIP (Table 4). There- 
fore, the detection of these O Vll features 
indicates the presence of a low ionization, 
high velocity component. According to our 
models (see §4.4 and Table 5) a component 
capable of reproducing the absorption lines 
in O VII but no absorption in O VIII or 
Ne X at 1365 km s^^, would need a low 
ionization parameter, but also a low col- 
umn density (log N// ~ 21) and therefore. 



would only produce observable absorption 
features by Ne VI- VII between 13.5 and 14 
A(ls — 2p transitions mainly), and N VI be- 
tween 23.1 and 23.8 A and at 24.9 A(ls-2p 
and Is — 3p transitions). Since the Ne VI- 
VII lines are weak and lie in a region dom- 
inated by the Fe L-shell transitions, they 
are undetectable. The N VI lines are also 
undetectable due to low signal to noise ra- 
tio and calibration uncertainties at long 
wavelengths. Therefore, the O Vll hues 
are really the only features observable from 
this component. This component would 
also produce absorption by O VI, and an 
absorber at roughly this outflow velocity 
has been already clearly observed in the 
UV (see §4.4). This UV detection gives 
us greater confidence in the reality of the 
high velocity LIP absorber, which is other- 
wise of borderline significance in the X-ray 
data. 

4.4. UV Absorbers 

Kraemer, Crenshaw & Gabcl (2001) and 
Gabcl et al. (2003) have presented the 
UV spectrum of NGC 3783. Absorption 
hues due to Civ, Nv, and Ovi are clearly 
present in the HST and FUSE spectra, 
with each line showing four different com- 
ponents corresponding to four outfiow ve- 
locities. One of these components (out- 
fiow 1027 km s~^) is too weak to be 
detected with X-ray instruments (see Fig. 
11). However, the other three components 
are consistent with the outfiow velocities 
of the X-ray LIP O Vll absorbers (at 750 
km s'^ and the 4 O Vll lines at 1345 km 
s~^). Furthermore, the ionization degree 
of the LIP X-ray components can also pro- 
duce UV features (in contrast with the 
HIP, which is almost transparent in the 
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Fig. 10. — Low ionization phase decom- 
posed in three subcomponents matching 
the FWHM, outflow velocity, and H col- 
umn density of the absorbers predicted in 
the UV band (Kraemer et al. 2001). The 
ionization parameters, however, are differ- 
ent (see §4.4). The model is presented only 
in the 15.5-19 A range. In other ranges, 
there is no significant difference with the 
model presented in Fig. 3. As can be ob- 
served, the two observed components of O 
VII are predicted. The bottom panels show 
the 3 most intense O Vll lines in velocity 
space. The model has been shifted up 10 
counts for clarity.^ 

UV band due to its high ionization degree). 
Motivated by this, and with the aim of ex- 
ploring the connection between UV and X- 
ray absorbers, we decided to model the X- 
ray data for comparison with the models 
by Kraemer et al. (2001). 

To begin with, we considered model 
parameters exactly as in Kraemer et al. 
(2001): the total equivalent column den- 
sity Nff, ionization parameter U, three dif- 
ferent velocity components for each line, 
with FWHM and outflow velocity Vout as 
for the UV lines (see Table 5). We found 



that the predicted strengths of X-ray lines 
do not match the data (see Fig. 12a). 
Then we changed only one input parame- 
ter, the ionization parameter U. We fitted 
the X-ray spectrum, leaving U free pa- 
rameter, but with all other parameters as 
given above. A good fit to the data was ob- 
tained, as shown in figure 10, for the 15.5 
to 19.0 A range (in other ranges, there is 
no significant difference with respect to the 
model presented in Figure 3, as discussed 
in §4.3.2). However, the best fit value of U 
is smaller by at least a factor of 3, than that 
found by Kraemer et al. for the UV spec- 
trum (Table 5). In principle, this can be an 
effect of the different input continuum that 
we used, so we repeated the analysis by us- 
ing exactly the same SED as in Kraemer et 
al. (i.e. with a value of F = 1.8) above 0.6 
keV and no high/low energy cutoffs). The 
resulting U still did not match the Krae- 
mer et al. value. The degree of ionization 
found by Kraemer et al. is too high to fit 
the X-ray absorber. A naive interpretation 
of these facts could be that the UV and X- 
ray absorbers are two different phenomena, 
independent of each other. This, however, 
is not true. As discussed below in §5, the 
high values of U inferred from the UV data 
alone, would produce strong absorption by 
the UTA in the region around 16 A and 
this is not consistent with the observations. 
Figure 12a illustrates this point clearly. 

The excellent level of agreement be- 
tween the data and our new model with 
three different velocity components for 
each line (as in the UV) is noteworthy. 
This strongly indicates a common origin 
of UV and X-ray absorbers, as we have 
been suggesting for years (Mathur et al. 
1994, 1995, 1997, 1998, 1999; Monier et al. 
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Table 4 

Most intense unblended absorption lines 



Equivalent Width (mA) 
Ion Wavelength (A) HIP LIP Total Measured 



Sixiv 


6.182 


16.1 


< 0.5 


16.1 


17.7 


Sixiii 


6.648 


20.4 


< 0.5 


20.4 


17.6 


Mgxii 


7.106 


7.3 


< 0.5 


7.3 


7.3 


Alxiii 


7.173 


3.7 


< 0.5 


3.7 


4.9 


Alxii 


7.757 


3.8 


< 0.5 


3.8 


3.4 


Fexxiii 


8.304 


5.1 


< 0.5 


5.1 


11. 7^= 


Mgxii 


8.421 


24.7 


< 0.5 


24.7 


27.3 


Fexxii 


8.714 


4.5 


< 0.5 


4.5 


5.2 


Nex 


10.239 


20.7 


< 0.5 


20.7 


23.0^ 


Fexix 


10.816 


18.2 


< 0.5 


18.2 


18.6 


Fexxii 


11.770 


28.4 


< 0.5 


28.4 


41.0^= 


Nex 


12.134 


45.4 


< 0.5 


45.4 


a 


Fexxi 


12.284 


39.7 


< 0.5 


39.7 


ATA" 


Fexx 


12.576 


35.7 


< 0.5 


35.7 


34.7 


Fexx 


12.754 


9.2 


< 0.5 


9.2 


12.3 


Fexx 


12.824 


26.1 


< 0.5 


26.1 


a 


Fexx 


12.846 


38.9 


< 0.5 


38.9 


a 


Fexx 


12.864 


37.0 


< 0.5 


37.0 


a 


Fexix 


13.423 


26.1 


< 0.5 


26.1 


a 


Neix 


13.447 


31.5 


11.0 


42.6 


a 


Fexix 


13.462 


30.0 


< 0.5 


30.0 


a 


Fexix 


13.497 


36.4 


< 0.5 


36.4 


a 


Fexix 


13.518 


48.3 


< 0.5 


48.3 


a 


Fexix 


13.795 


28.5 


< 0.5 


28.5 


a 


Nevii 


13.814 


< 0.5 


42.3 


42.3 


a 


Nevi 


14.020 


< 0.5 


41.1 


41.1 


a 


NeVI 


14.047 


< 0.5 


51.4 


51.4 


a 


Fexviii 


14.208 


49.2 


< 0.5 


49.2 


a 


Fexviii 


14.208 


41.4 


< 0.5 


41.4 


a 


Fexviii 


14.373 


30.6 


< 0.5 


30.6 


a 


Oviii 


14.821 


22.3 


3.1 


25.4 


36.5 


Fexix 


14.961 


9.8 


< 0.5 


9.8 


9.5^ 


Fexvil 


15.014 


54.6 


< 0.5 


54.6 


33.7<= 



+1.9 
-1.9 
+2.9 
-2.7 
+1.9 
-2.8 
+2.1 
-2.0 
+3.1 
-2.5 
+4.3 
-4.1 
+3.8 
-3.4 
+3.7 
-2.3 
+3.8 
-3.2 
+4.8 
-3.9 
+8.8 
-7.4 



+8.3 

-7.1 

+10.2 

-8.7 

+3.8 

-4.3 



+15.1 

-12.6 

+6.4 

-5.6 

+17.6 

-15.2 
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Table 4 — Continued 



Equivalent Width (mA) 



Ion 


Wavelength (A) 


HIP 


LIP 


Total 


Measured 


Oviii 


15.176 


34.6 


6.5 


41.1 


38.7 


+14.4 
-16.2 


Fexvii 


15.261 


24.5 


< 0.5 


24.5 


19.8 


+6.5 
-5.9 


Oviii 


16.006 


51.0 


16.7 


67.7 


59.8^ 


+10.8 
-17.4 


Ovii 


17.086 


< 0.5 


24.5 


24.5 


b 




Ovii 


17.200 


< 0.5 


33.0 


33.0 


28.9^ 


+6.8 
-7.6 


Ovii 


17.396 


< 0.5 


43.4 


43.4 


38.9^ 


+7.8 
-11.2 


Ovii 


17.768 


5.0 


55.2 


60.2 


40.1^ 


+16.0 
-15.2 


Ovii 


18.627 


13.2 


71.0 


84.2 


46.3'' 


+16.9 
-14.2 


Oviii 


18.969 


89.0 


54.1 


143.1 


55.0'' 


+23.8 
-15.4 


Nvii 


20.910 


19.3 


12.8 


32.1 


29.5'' 


+23.3 
-21.4 


Ovii 


21.602 


49.5 


128.3 


177.8 


34.2'' 


+30.8 
-33.1 


Nvi 


23.277 


< 0.5 


24.8 


24.8 


a 




Nvi 


23.771 


< 0.5 


41.9 


41.9 


23.9 


+21.3 
-22.5 


Caxiv 


24.114 


25.2 


< 0.5 


25.2 


26.8 


+20.2 
-19.9 


Nvii 


24.781 


67.0 


54.9 


121.9 


66.5'' 


+37.2 

-34.4 


Nvi 


24.898 


< 0.5 


68.7 


68.7 


33.9'' 


+31.7 
-27.2 



^Not measured due to heavy blending with several absorption 
lines 

^EWs do not take into account blending with the correspond- 
ing emission line 

'^EWs do not take into account blending with other absorption 
line 
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2001, although with a one phase model; 
only high resolution spectra allows the 
identification of two components). The 
above exercise also highlights the danger 
of modelhng the absorbers using UV data 
alone. The UV value of the ionization 
parameter underpredicts the ionic column 
densities by as much as factors of ten. The 
most likely explanation of this discrepancy 
is that the UV lines are saturated and the 
saturation effects are often estimated in- 
adequately. This is illustrated beautifully 
in the case of NGC 5548 by Arav et al. 
(2003). They find that the Ovi column 
density measured from UV data by Broth- 
erton et al. (2002) is lower than the ro- 
bust lower limit obtained from the XMM- 
Newton data, by at least a factor of seven. 
When the effects of saturation are dealt 
with in a more sophisticated way by Arav 
et al. a good agreement is found. 

Saturation was the main point in the 
models of Mathur et al. who parameter- 
ized its effect using the simple technique 
of curve-of-growth analysis. The key point 
that we stress again here is that the X- 
ray data imply a large amount of UV ab- 
sorption, which cannot be hidden by claim- 
ing that the UV and X-ray absorbers are 
different. We will discuss the relationship 
between the UV and X-ray absorption in 
NGC 3783 in a forthcoming paper. 

5. Evaluation of the Spectral Fit- 
ting Results 

The self-consistent absorption model we 
present in this paper is the first to be car- 
ried out for the 900 ksec spectrum of NGC 
3783. The extraordinary signal to noise 
ratio of the data has allowed us to obtain 
important constraints that make the mod- 



elling results more reliable. This is an im- 
portant point since without the restrictions 
imposed by the data (see below), small 
changes in the fitting can give rise to com- 
pletely different physical interpretations of 
the medium (see §5.5). 

Our model presents a surprisingly sim- 
ple picture: two kinematically indistin- 
guishable gas outflows with quite different 




0.1 1 10 



u 

Fig. 11. — Equivalent hydrogen column 
density (^h) vs. ionization parameter 
{U = riphot/ne), showing observational 
lower limits for two UV lines (to a 3cr 
level): CIV (HST/STIS), OVI (FUSE), 
and two X-ray lines: OVII, OVIII {Chan- 
dra/LETG). The limits are calculated for 
lines of width, 6=100 km s~^, a UV flux 
f(1032A)=5xlO-^%g cm-^s-^ 20 ksec 
HST and 50 ksec FUSE exposures; an X- 
ray flux of f(0.5-2 keV)=2 mCrab, and a 
100 ksec LETG exposure. The shaded re- 
gion shows where the OVI measurement is 
unsaturated (i.e. lies on the linear part of 
the curve of growth).^ 
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Table 5 

Parameters found for U, assuming UV absorber values for Nh, FWHM, 

AND Yout- 



Parameter 




Comp 1 


Comp 2 


Comp 3 


Log Nh {cm-^Y 


fix 


21.08 


20.81 


21.18 


FWHM (km s"^)" 


fix 


193 


170 


280 


Yout (km s-i)" 


fix 


1365 


548 


724 


Log U(UV)" 




-0.11 


-0.1 


-0.19 


Log U (X-ray)'' 


free 


-1.20±0.09 


-0.69±0.21 


-0.99±0.11 



'^Values predicted in tlie UV by Kraemer et al. 2001. 

*^Thc only free parameter in tliis model. Calculated assuming 
N^^, FWHM, and Vout values from the UV. 



photoionization equilibrium temperatures 
and two widely separated ionization pa- 
rameter values. The model also suggests 
that no third component with high H col- 
umn density and intermediate U is consis- 
tent with the data (see below) . Six free pa- 
rameters are sufficient to model > 100 ab- 
sorption features and to correctly predict 
the EW of the 23 most intense lines that 
have little or no blending. Both compo- 
nents are consistent with solar abundances. 
Our fits are sensitive to departures of the 
Mg/0, Fe/0, Fe/Mg, Ne/0 ratios by any 
more than a factor of ~2 from solar, since 
transitions from these elements are in re- 
gions of high S/N. These data are not sen- 
sitive, however, to the C/N relative abun- 
dance. A component with low U and Nh at 
1345 km s~^ is needed to model an asym- 
metry in 4 O VII lines. A fourth compo- 
nent is evidenced by the presence of two 
Fe XXV absorption lines at 1.57 and 1.85 



A (identified by Kaspi et al. 2002). How- 
ever, such component is even more highly 
ionized, and only contributes with a few 
lines to the overall spectrum. Almost all 
these features are below 4 A, and are only 
detected by the HEG configuration of the 
HETG (not modelled here). This compo- 
nent is not well constrained (particularly 
by the MEG band), and we defer this anal- 
ysis to a forthcoming paper. 

Absorption by iron and oxygen is partic- 
ularly important in the spectrum of NGC 
3783. Even though some of their fea- 
tures arc blended, these elements can be 
regarded as the main "tracers" of the gas, 
and therefore, can be used to evaluate the 
validity of a given model. In our case, an 
important result is that while O Vlll and 
Fe xvil-xxil arc the main features of the 
hotter phase (sec Fig. 7), the signature of 
the cooler phase can be attributed to the 
UTA and the O Vll hues (see Fig. 8). Only 
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2% of the O VII, but 90% of the O viii, 
arises in the HIP (Ni/p(0 vil)^ 1.5 x 10^^ 
cm"^ and Nl/p(0 VIII)?« 1.2 x 10^'' cm~^ 
vs. Nhip{0 VIi)^ 3.4 X 10^^ cm-^ and 
Nhip{0 vm)^ 1.2 X 10^^ cm-^). This 
is noteworthy, since traditionally, the O 
VII absorption has been associated with 
that of O VIII and warm absorbers (e.g. 
George et al. 1998 and references therein). 
Otani et al. (1996) have previously sug- 
gested two different sources of O Vll and 
O VIII (for MCG-6-30-15), but their phys- 
ical interpretation is completely different 
than the one presented here, with the two 
absorption components at completely dif- 
ferent regions. In our model, O Vll is 
consistent with a cooler medium produc- 
ing significant amounts of O VI and other 
lower ionization species. Without the con- 
straints introduced by the UTA this dis- 
tinction could not be realized. We also 
note that the O column densities produce 
significant K-cdgc absorption (r^ = 0.33 
for O VII and ta = 0.13 for O vill). 

Our model predicts more than 265 ab- 
sorption lines with EWs larger than 1 mA. 
As can be appreciated in Figure 3, most 
of them fit the data in a quite satisfac- 
tory way. Table 4 and Figure 9 also show 
a good agreement between the measured 
and predicted EWs for unblended lines. 
The highly ionized phase produces nearly 
all the features between 15.3 A and 8.5 
A (Fig. 7), where the spectrum becomes 
dominated by the Fe L-shell lines, plus the 
contribution of high ionization lines arising 
from other ions, such as O, N, and Mg. The 
predicted form of the UTA by the cooler 
phase also shows a good agreement with 
the data. We notice, however, that individ- 
ual fines at around 16.15 A and 16.35 A are 



not predicted by the model. These features 
are part of the UTA (as can be observed in 
Figure 4 of Behar, Sako, & Kahn; 2001), 
but the abbreviated data cannot reproduce 
them. Our model does not yet include data 
for the O VI inner shell transitions. How- 
ever, we notice that these features can be 
observed in the spectrum, particularly the 
most prominent one at 22.04 A. Unfortu- 
nately, the signal to noise ratio at these 
wavelengths is too low to confirm their de- 
tection, even at the 1 a level (see also 
Fig. 11). This is also the case of the fea- 
ture at 22.33 A, which could arise from O 
V. However, the presence of these features 
in the spectrum of NGC 3783 has been 
established through observations with the 
XMM-Newton (Blustin et al. 2000). This 
is an important point, since the LIP in our 
model predicts that a substantial fraction 
of the O will exist in these lower charge 
states (N(0 Vl)^ 7 x 10^^ cm'^ and N(0 
Y)Ri 2 X 10^^ cm-2). 

There is another region where a discrep- 
ancy deserves comment. Between 6.75 and 
7.0 A, two lines clearly show larger EWs 
than the predictions of the model. These 
intriguing lines may be identified provi- 
sionally as Si x(A6.850) and Si xi(A6.775). 
To reproduce them, either the high ioniza- 
tion phase would have to be much cooler 
(T < 2 X 10^ K) with a much smaller value 
of the ionization parameter {\ogU < 0.5), 
or the low ionization phase would have to 
be much hotter (T > 9 x 10^ K) with a 
much larger value of the ionization param- 
eter (logC/ > —0.2). As will be discussed 
in §5.1, such values are not consistent with 
other features observed in the spectra. The 
presence of a third absorber with interme- 
diate U is also unhkely (see Fig. 12d). 
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A possible explanation for the difference 
is that our atomic database is incomplete, 
containing data for only the most intense 
lines of Si (lines with oscillator strengths 
larger than 0.1), so that the observed fea- 
tures are in fact blends of several Si x and 
Si XI lines. Alternatively, the atomic data 
used to calculate the ionization balance for 
iron and silicon might be inconsistent. 

Our model uses the abbreviated data set 
to fit the 15-18 A UTA. Since this is the 
main feature of the cooler phase, it is im- 
portant to know the effects that this lim- 
ited data could have on our results. We 
believe the approximation is sufficiently re- 
liable in this case for two reasons. First, 
we are using turbulent velocities > 100 
km s~^ and for all ions of iron the column 
density is lower than 10^^ cm~^ (which ac- 
cording to Behar, Sako, & Kahn (2001), 
are the limits within which the approxima- 
tion is valid, see §3). Second, our results 
are in qualitative agreement with those by 
Blustin et al. (2000) who modelled the 
UTA with the complete data set, and find 
a cooler phase consistent with ours (see 
§5.5). Therefore, it is unlikely that an over- 
estimate of the absorption is present, but 
this still has to be confirmed. A quanti- 
tative estimation of the uncertainties in- 
troduced by the abbreviated data will be 
presented in a forthcoming paper. 

5.1. Ionization constraints on two, 
three, and multi-component 
models 

Since it is a wide, well resolved feature, 
the Fe M-shell UTA sets tight restrictions 
on the ionization level of the gas (see Fig. 
12). As can be observed in Figure 13, dif- 
ferent values of the ionization parameter 
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Fig. 12. — Panel (a): Low ionization component 
decomposed in three subcomponents matching all 
the parameters predicted for the UV absorbers 
(Kraemer et al. 2001). This fit shows several dis- 
crepancies with the data. Middle panels: Differ- 
ent values for the ionization parameters for our 
two phase absorber showing inconsistencies with 
the observed UTA. Panel (b): logU=0.5 for the 
high ionization phase. Note that the absorption 
from Si XIV is underestimated. Panel (c): In this 
case, the ionization parameter of the low ioniza- 
tion component has been set to logU=-0.2. Panel 
(d): A third absorber has been included in our 
model with logU=-0.2 and logN H=21.7. The in- 
clusion of a new component could account for the 
missing Si lines, but would overpredict the absorp- 
tion by the UTA, therefore a flow with continuous 
variation of U is ruled out.^ 



lead to unambiguously different shapes for 
the UTA. For instance, logU = —1.8 will 
produce a feature from Fe lll-vil between 
16.5 and 17.5 A, with the peak of absorp- 
tion around 17 A (Fig. 13a). Lower val- 
ues of logU will give rise to a narrower 
UTA, centered at wavelengths larger than 
17 A. On the other hand, larger values 
for the ionization parameter will shift the 
peak of absorption to shorter wavelengths. 
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Fig. 13. — UTA model predictions for dif- 
ferent values of the ionization parameter. 
The shape of the UTA is very well defined 
by the degree of ionization, which results in 
important constraints to the physical con- 
ditions of the ionized gas (Behar, Sako & 
Kahn, 2001). 

and will produce a wider shape. A log- 
arithm of the ionization parameter equal 
to -0.78, hke the one predicted from the 
LIP in OTU^ model, produces absorption be- 
tween 15.7 and 16.9 A, with an absorption 
peak at around 16.4 A (Fig. 13b). A fur- 
ther increase in the ionization parameter 
value will further shift the peak of absorp- 
tion toward smaller wavelengths and will 
widen even more the shape of the UTA 
(Fig. 13c). 

The UTA can thus be used to estimate 
a lower limit to the ionization parameter 
of additional ionization states in the ab- 
sorber. The spectrum of NGC 3783 shows 
no trace of the inner shell Fe XIV-XVI M- 
shell transitions, and hence most of the 



iron in the HIP has to be at least 15 times 
ionized. The fractions of Fe XV and Fe XVI 
for this component are less than ~ 0.003 
and ~ 0.001 respectively. 

The ionization parameter of the HIP is 
not consistent with a smaller value capable 
of reproducing the Si absorption features 
between 6.75 and 7.0 A, since it would pro- 
duce a significant fraction of Fe XII-XV, 
and a UTA similar to the one presented 
in the lower panel of Figure 13, a feature 
not observed in the spectrum (see section 
§5.1). Furthermore, a larger value would 
undcrpredict the absorption by Si XIV and 
other high ionization species. The result- 
ing inconsistencies can be observed in Fig- 
ure 12b. The low ionization phase is sim- 
ilarly well constrained. Figure 12c shows 
the effects of a larger ionization parame- 
ter. In this case, the absorption by the 
UTA at wavelengths larger than 16.5 A is 
underpredicted by > 50 %. 

The inclusion of another absorber with 
intermediate ionization parameter and 
comparable H column density is also in- 
consistent with the data. A third compo- 
nent, capable of reproducing the Si lines, 
would overpredict the UTA absorption be- 
tween 15.3 and 16 A by > 75%. Such an 
absorber would need an ionization param- 
eter logU ~ —0.2 and a column density 
logNij 21.7. Figure 12d presents the ef- 
fects of including a third absorber. These 
constraints strongly suggest that no in- 
termediate values of U with high column 
density are present. However, the pres- 
ence of the Si X and xi absorption features 
not accounted for in our model docs not 
allow us to effectively rule out this possi- 
bility. Nevertheless, the use of better data 
for both the UTA and the Si inner shell 
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transitions may account for the discrepan- 
cies. We cannot rule out the presence of 
additional low column density components 
or a bumpy continuous distribution with a 
small contribution to the absorption by low 
column density gas. However, such puta- 
tive gas would have to have small column 
densities (logN// < 21 [cm~^]). Hence, 
from the above analysis it is clear that 
only two dominant absorbing components 
are required by the data in the MEG wave- 
length range. 

As a final remark, we would like to stress 
that the discussion presented in this sec- 
tion is not dependent on the fact that we 
are using the abbreviated data to model 
the UTA. This approximation affects only 
sharp narrow edges. Figure 4 of Behar, 
Sako & Kahn (2001) presents three UTA 
models produced by media with different 
ionization parameters carried out with the 
full set of Fe M-shell transitions. A close 
inspection of their figure leads to the same 
conclusions presented here. 

The physical conditions derived for the 
ionized absorber in NGC 3783 are thus well 
constrained. 

5.2. Line Saturation 

Our model predicts saturation of sev- 
eral absorption lines (see for instance Fig. 
6). This saturation is not evident to the 
eye due either to heavy line blending or 
to emission and absorption line filling (see 
§5.3). Figure 14 shows a comparison of the 
EWs predicted from our model calculated 
with a gaussian and a Voigt profile. Fewer 
than 20 lines show saturation. However, 
as shown in the figure, the most intense 
absorption lines from the LIP are satu- 
rated. From Figure 8, it is evident that the 
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Fig. 14. — Model predicted EWs calcu- 
lated using a Voigt profile vs. calculated 
using a gaussian profile. The absorption 
lines from the HIP and the LIP are pre- 
sented separately. Deviation from the 1:1 
identity line indicates strong saturation.^ 

LIP contributes with fewer absorption lines 
(but is responsible for most of the photo- 
electric absorption) than the HIP. So, al- 
though present in only a few absorption 
lines (the most intense ones) saturation be- 
comes a significant effect. 

5.3. Absorption and Emission Blend- 
ing 

As can be observed in Figure 3, the 
model predictions are in fairly good agree- 
ment with the data for the most intense 
O VII and O vill absorption and emission 
lines. However, an inspection of these lines 
from Table 4 shows that the EWs derived 
by the model have systematically larger 
values than those measured from the ob- 
servations. We believe this effect is a re- 
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suit of the filling of the absorption fines by 

the emission features. In Figure 6, we illus- 
trate an example of the blending of absorp- 
tion and emission lines. Therefore, one re- 
ally needs to take into account absorption 
and emission together to correctly infer the 
physical conditions of both the absorber 
and the emitter. In our case, the emission 
was not inferred in a self-consistent way, 
although, it was constrained by the ab- 
sorption. Therefore our emission fits (Ta- 
ble 2) can be used as a rough estimate of 
the emission line properties of NGC 3783. 
An inevitable conclusion is that without 
modelling of absorption and emission line 
blending, measurements of the observed 
EWs from the data will underestimate the 
actual values. This effect can be observed 
in Figure 15a, where we present a com- 
parison of the emission line EWs with and 
without the contribution from absorption. 
This figure also shows that a factor of 2 
for the filling of the absorption lines by 
the emission lines is present in almost all 
the resonant transitions (the few forbidden 
transitions show no effect, as expected). 

Additional complexity is introduced to 
the absorption/emission blending problem 
by the continuum level: if this level is not 
modelled globally, inconsistent results can 
be obtained. In Figure 15b, we compare 
the emission properties from this paper 
with those deduced by Kaspi et al. (2002), 
where no filling correction was applied, and 
only a local continuum was fitted. For 
truly unblended lines (such as O Vll forbid- 
den and intercombination lines) , the values 
reported by Kaspi et al. are larger than 
ours, indicating a lower continuum level. 
By contrast, some emission lines blended 
with an absorption component have lower 
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Fig. 15. — EW comparison for emission lines. 
Crosses within a circle stand for O vn forbid- 
den and intercombination lines, where no ab- 
sorption and emission filling is expected. Up- 
per panel: EWs comparison between measurement 
taking and without taking into account blending 
and absorption/emission filling model predictions. 
Lower panel: EWs from Kaspi et al. (2002) vs. 
EWs from this paper. In this case, since the 
emission/absorption blending is not taken into ac- 
count, and only a local continuum is fitted, incon- 
sistent results arc obtained. The solid line has a 
slope equal to 1. The dashed lines represent devi- 
ation factors of 1/2 and 2. 

EWs in their study (in some cases by a 
factor 2), a signature of absorption filling. 
Even more, the fitting of a local contin- 
uum can yield highly misleading results. 
For instance, Kaspi et al. estimated the 
EW of the resonant Ne lx(A13.699) line as 
78.1 ±12. 6 mA, while we computed a value 
of only 15.3lg'g^ mA. Although this line is 
filled by absorption, Kaspi et al. found 
a much larger EW than we do, because 
they assumed that most of the spectral fea- 
tures near this wavelength were produced 
by emission (i.e. they fitted a low contin- 
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uum level). As can be seen in Figure 3 
the continuum level at this point from our 
global fit is rather high, and most of the 
structures there are produced by absorp- 
tion, with only a small contribution from 
this emission line. This is a common prob- 
lem also faced in analyses of the Lyman a 
forest (e.g. Bechtold et al. 2002). 

5.4. Spectral Energy Distribution 

Mathur, ct al. (1994) have shown 
the importance of including the particu- 
lar SED of an AGN, rather than a typical 
SED. In the latter case, incompatible phys- 
ical conditions may arise in the description 
of a photoionized cloud. As discussed be- 
fore (§4.1), the SED we used to calculate 
the ionization balance of the gas is only 
an approximation to the actual (unknown) 
SED. This approximation is the best that 
can be obtained currently for the spectral 
shape of NGC 3783. This point is impor- 
tant, since the values obtained for the ion- 
ization parameters are dependent on the 
continuum shape, and moreover, so is the 
difference between the ionization parame- 
ter of the two phases. 

Our chosen SED (Fig. 5) produces 
rapid changes of the ionization balance 
with small changes of the ionization pa- 
rameter, and therefore the ratio of U be- 
tween the two phases is small when com- 
pared to other SEDs (for instance the SED 
used by Kaspi et al. 2001). To ilhistrate 
this point, we have recalculated the ioniza- 
tion parameters using the Mathews & Fcr- 
land (1987) continuum (shown in Fig 5). 
The values obtained for this SED for the 
low and high ionization component were 
logf/ = —0.05 and logf/ = 1.81. These 
values are larger than those predicted with 



our better chosen SED, and yield a fac- 
tor K, 72 between the ionization parameter 
of the phases, twice as large as we obtain. 
Therefore, any comparison between models 
should account for the effects of the differ- 
ent SEDs used. In our case, a factor of 
only 35 in the ionization parameters pro- 
duces dramatically different ionization de- 
grees, with a factor ^ 37 between the tem- 
perature of the two phases. However, using 
the Mathews & Ferland continuum gives 
T(HlP)=8.1xlO^ K and T(LlP)=2.1xlO^ 
K, and so a difference factor in tempera- 
ture of ~ 40, all of which are close to the 
factor found with our SED. 

5.5. Previous Models 

Four different studies have been pub- 
lished in the past for NGC 3783: (1) Kaspi 
et al. (2001) modelled a much lower qual- 
ity spectrum, the first HETG- C/iandra ob- 
servation of NGG 3783; (2) Blustin et al. 
(2002) modelled also lower quality data, 
a 40 ksec exposure taken with the XMM- 
Newton; (3) de Rosa et al. (2002) modelled 
low spectral resolution BeppoSAX obser- 
vations; and (4) Kaspi et al. (2002) ana- 
lyzed the same 900 ksec HETG spectrum 
with an empirical approach. In the follow- 
ing, we review these analyses and contrast 
them with our results. Table 6 summarizes 
the results found with the different models 
(including the one presented in this paper) . 

5.5.1. The model by Kaspi et al. 2001 

These authors established for the first 
time that two absorption systems were re- 
quired to fit the wide range of ionization 
species found in the data. They estimated 
the continuum SED, including (1) spec- 
tral curvature resulting from bound-free 
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Table 6 

Comparison with previous models 





r 


Norm." 


log(U^f) 


log(UgF) 


log(N^^^) 


log(Ng^^) 


z(x)/zo 


This paper 


1.53'' 


.011 


-2.77 


-1.23 


21.6 


22.2 


1 


Kaspi et al. 


1.77 


0.015 


-1.75 


0.75 


22.2 


22.2 


1 


Blustin et al 


1.53 


.018 


0.3^ 


2.4:" 


20.7 


22.45 


lO(Fe) 


de Rosa et al. 


1.83'^ 


0.019 




0.4 




22.3 


1 



^In photons keV ^ cm ^ s ^ at 1 keV. 

'^The continuum includes the contribution from a black body component with kT=0.1 
keV, and norm=0.0002 L^g/Df^, where L39 is the source luminosity in units of 10^^ erg s~^ 
and DiQ is the distance to the source in units of 10 kpc. 

^The value reported is log^, where ^ = L/r^n, L the luminosity of the source, r the 
distance from the source to the cloud, and n the cloud density. A conversion to Uqx was 
not possible since the authors did not provide the SED. 

'^The continuum includes the contribution from a black body component with kT=0.21 
keV, and norm=0.00045 L^g/Dl^. 
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absorptions due to a single component 
ionized absorber and (2) "Compton re- 
flected" continuum. Using this SED, they 
then modelled the absorption assuming the 
same turbulent and outflow velocities for 
the two components (300 km s~^ and 610 
km s~^ respectively). They found a photon 
index F = 1.77 and a normalization value 
A(lkeV)=0.0151 photons keV"^ cm'^ g-^ 
While for their low ionization compo- 
nent they found logUoxYGEN — —1.745, 
for the high ionization one they obtained 
logUoxYGEN = —0.745. In both compo- 
nents log Nj^[cm~^] was estimated to be 
22.2. 

A comparison between this model and 
ours (see Table 6) indicates a difference of 
0.24 between the photon indexes deduced 
(r = 1.53 in our case). Compton scat- 
tering (not included in our model) is un- 
likely to contribute a lot to this difference; 
Kaspi et al. (2002) estimated it to be 
< 15%. Rather, the source of the discrep- 
ancy most likely arises from Kaspi et al. 
(2001) basing their column density mea- 
surements only on the photoelectric edges. 
The 900 ksec spectrum clearly shows that 
the two deepest edges, O VII and O VIII, 
are masked by Fe transitions: O VII is hid- 
den in the Fe M-shell UTA (not modelled 
by Kaspi et al.), and O Vlll is hidden in 
the Fe L-shell lines. Since the contribution 
of these lines was not considered, the edge 
depths were overpredicted (see Fig. 4), 
which, in turn, yielded an overprediction 
of the photon index. This effect has reper- 
cussions on any measurements based solely 
on low resolution CCD spectra. In the 
case of NGC 3783, based on ASCA data, 
Reynolds (1997) and George et al. (1998) 
found Tx = 1.2 for O Vll and ta = 1.4 



for O VIII, and thus a factor at least 4 

times larger than the edges predicted by 
our model. Figure 4 shows a compari- 
son between models including only photo- 
electric absorption and models including 
both photoelectric and resonant absorp- 
tion. Therefore the depths and the column 
densities calculated through low resolution 
spectra, should be considered only upper 
limits to the actual values. In objects hke 
NGC 3783, where a deep UTA is present, 
the ovcrcstimation is likely to be signifi- 
cant. We do not yet know the importance 
of this effect in other objects. 

To make a comparison of the ionization 
parameters, we calculated logUoxYGEN 
(the ionization parameter in the range 
0.538 and 10 keV, George et al. 2000) 
for the continuum we used. The values we 
obtained were -2.77 for the LIP and -1.23 
for the HIP, factors of 10 and 3 lower than 
Kaspi et al. (2001) values. While our HIP 
has an ionization degree between the high 
and the low ionized components of Kaspi 
et al., our LIP is much less ionized than 
either. In both models, the low ionization 
phase is responsible for the O Vll absorp- 
tion, but the predicted physical conditions 
of the gas arc completely different. The 
difference in ionization degree is explained 
by the omission of the UTA in their model. 
As discussed in §5.1, the ionization param- 
eter of their low ionization component is 
too high to reproduce this feature. 

5.5.2. The model by Blustin et al. 2002 

As in the previous case, Blustin et al. 
(2002) needed two absorption components 
to fit the XMM-RGS data. Blustin, et 
al. modelled the continuum with a simple 
power law with F = 1.54 and a normal- 
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ization value A(lkeV)=0.0175. The nor- 
malization is larger than what we find, be- 
cause they have not included a blackbody 
component in their model. The photon in- 
dex, however, is consistent with the one 
we found (see Table 6). Their strategy for 
modelling the absorption was to measure 
first the contribution of each ion with the 
configuration slab of the code SPEX, and 
then to model the data for a two phase ion- 
ized absorber using the configuration xabs. 
The authors estimated the outfiow velocity 
as 800 km s~^, and fixed a velocity width 
(rms) of 300 km s~^. This width refers to 
the gaussian sigma, and corresponds to a 
Doppler velocity (turbulent plus thermal, 
see §3) fti 425 km s~^. The values predicted 
by their model were log^ = 2.4 erg cm s~^ 
and Nh = 2.8 x 10^^ cm~^ for their high 
ionization phase, and log^ = 0.3 erg cm 
s"^ and Nh = 5.4 x 10^° cm~^ for the low 
ionization one (where ^, defined as L/nr"^, 
is the ionization parameter). A direct com- 
parison between ionization parameters is 
not possible since the SEDs used are differ- 
ent (see §5.4). However, our results are not 
qualitatively in contradiction with theirs. 

Nevertheless, there are several quantita- 
tive differences that lead to different inter- 
pretations of the data: (1) While we were 
able to fit the data assuming solar abun- 
dances, their low ionization phase requires 
10 times the solar abundance of iron to re- 
produce the observed features. This value 
seems to be rather high. (2) The H column 
density of their cooler phase IS ^ 8 times 
lower than ours. (3) By contrast with our 
model, in theirs the hotter component has 
a significant contribution to the absorption 
produced by O Vll, due to the high column 
density of this component (and the low col- 



umn density of their cooler component) . 

These quantitative differences can be 
understood from the different continuum 
fits applied to the data, i.e. due to the in- 
clusion of a soft thermal component in our 
model. The presence of a soft thermal com- 
ponent requires deeper absorption at long 
wavelengths to match the data. There- 
fore, a cooler phase capable of accounting 
for this absorption needs to have a high 
H column density. However, this would 
also result in strong saturation in the C, 
N, and O lines. Since, as stated by Blustin 
et al., the abundance of Fc is effectively 
measured against the abundance of these 
ions, omitting the soft excess in the pres- 
ence of a deep UTA would require an over- 
abundance of Fe relative to other elements. 
Therefore, Blustin et al. concluded that 
iron was around one order of magnitude 
more abundant than oxygen (relative to so- 
lar values) and, accordingly, that the H col- 
umn density was around one order of mag- 
nitude smaller in the low ionization phase 
because, at the time, they did not find evi- 
dence for a soft excess. However, they em- 
phasize that if saturation is present (along 
with the presence of a soft excess), the iron 
abundance could be much smaller. This is 
in agreement with our results, which re- 
quire Fe/0 to be within a factor 2 of solar 
and strong saturation. The inclusion of a 
thermal component, therefore, allowed us 
to fit the data with solar abundances. 

The existence of a soft excess in the 
spectrum of NGC 3783 has been further 
evidenced by the X-ray variability proper- 
ties seen with new XMM data. The over- 
all 0.2-10 keV spectrum gets softer as the 
source gets brighter, consistent with a vari- 
able soft excess, as seen in other Seyfert 
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galaxies (A. Blustin, private communica- 
tion) . 

5.5.3. The Model by de Rosa et al. 2002 

BeppoSAX observed NGC 3783 5 days 
in June of 1998. De Rosa et al. 2002 re- 
ported this observation and modelled the 
data. Only modest flux variations (~ 
20%) were detected. They found a pho- 
ton index F = 1.83 and a normalization 
value A(lkeV)=0.019 photons keV^^ cm~^ 
s~^. A thermal component with kT— 0.21 
keV and normalization 4.5x10^^ L,39 D\q 
(where L39 is the source luminosity in units 
of 10^^ erg s^^ and Diq is the distance to 
the source in units of 10 kpc) was neces- 
sary to fit the spectrum. They also found 
evidence for the ionized absorber. The re- 
sults of their fitting are presented in Table 
6. As in the case of Kaspi et al. (2001), 
the differences between fittings can be ex- 
plained by the fact de Rosa et al. used 
only the absorption edges (due to the low 
resolution of their data), and not the lines, 
in their analysis. 

5.5.4. Analysis by Kaspi et al. 2002 

Kaspi et al. (2002) used a different 
analysis strategy and did not fit the data 
with a theoretical model. They empirically 
identified the lines and measured the EWs 
directly from the spectrum, using appar- 
ent line free zones to estimate the contin- 
uum level, but only locally. The appeal of 
this approach is that it produces directly 
measured quantities before getting into de- 
tailed modelling. The disadvantages are: a 
confiating of blends into single "features" , 
the ignoring of the mutual cancelling of 
emission and absorption, and an inability 
to derive a global picture of the physical 



conditions. 

There are some discrepancies between 
their line identifications and the predic- 
tions from our model. The most important 
ones are listed in Table 7. The differences 
can be attributed mostly to the inclusion in 
our calculations of Ni and Ca lines, as well 
as more than 1200 Fe transitions arising 
from the ground state. According to Kaspi 
et al. (2001) their models include several 
hundred. The authors tentatively identi- 
fied a 14±16 mA absorption line around 
21.45 A (in the rest frame system, Fig 3) as 
a local component of O VIl(A21.602). This 
feature would be produced by intergalactic 
absorption at zero redshift, or by Galac- 
tic absorption (see Nicastro et al. 2002). 
Our model indicates that this line is pro- 
duced by the high ionization phase, and 
corresponds to a Ca xvi(A21.45) transi- 
tion (the difference in wavelength is due 
to the outflow velocity of the absorber). 
Therefore, any putative contribution from 
a local medium cannot account for an EW 
larger than 9.3 mA (at la level), consis- 
tent with the results found for the local fll- 
ament by Nicastro et al. (2002, 2003). An 
EW comparison between Kaspi et al. and 
our measurements indicate a good agree- 
ment for all unblended absorption lines (14 
cases) and for some blended features (11 
out of 17 cases). The differences found 
can be ascribed to the different contin- 
uum levels. While we fltted the contin- 
uum derived from the global spectrum us- 
ing our model, they measured the contin- 
uum locally and in an empirical way. We 
notice that, as discussed in §5.3, several 
absorption-emission blended lines present 
systematic differences between measure- 
ments and our model predictions. 
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Table 7 

Comparison with Kaspi et al. (2002) 



A (A) Ion Name and Rest Frame A (A) identification 



Observed 


Kaspi et al. 


This paper 


10.126 


Fexvii(10.112) 


Fexvil(10.112), Fexix(10.119), Fexvil(10.120) 


10.524 


Fexvil(10.504) 


Fexvii(10.504), Fexviii(10.537) 


12.436 


No identification 


Nixix(12.435) 


12.560 


Fexx(12.576) 


Fexx(12.526) 


12.592 


Fexx( 12.588) 


Fexx(12.576,12.588). Fcxix(12.538) 


13.612 




Fexix(13.645,13.643) 


13.822 


Fexix(13.795) 


NeVll(13.814),Fexlx(13.795) 


14.269 


Fexviii(14.256) 


Nev(14.239),Fexvill(14.256) 


15.584 




Fexviii(15.625) 


21.466 


Ovil(21.602) zero redshift 


Caxvi(21.450) 


23.783 




Nevi(23.771) 


24.124 




Caxiv(24.114) 



6. Absorber Physics Discussion 

6.1. Two phases of the same medium? 

The simple nature of our two compo- 
nent model, with indistinguishable kine- 
matics and very different ionization pa- 
rameters, suggests that the absorption ob- 
served in the NGC 3783 spectrum may 
arise from two phases of the same medium. 

First, two separate components seem to 
be the best way of describing the absorbing 
gas, since the discussion presented in §5.1 
and the plots in Fig. 12d strongly suggest 
that a continuous medium, i.e, a medium 
with a smooth distribution of temperatures 
from the HIP to the LIP (as proposed by 
Krolik & Kriss, 2001) is unlikely to pro- 
duce the well defined UTA absorption fea- 



ture observed. Rather, such a medium 
would produce less deep, continuous, ab- 
sorption from s:il5to~17A.A continu- 
ous medium also requires strong iron L and 
K edges (Krolik & Kriss 1995), features not 
detected in the spectrum. 

The presence of a two phase medium 
(and not a continuous medium) is fur- 
ther evidenced by the fact that the two 
components are in pressure equilibrium: 
the gas pressure (P) of each phase is pro- 
portional to the electron density (ue) and 
the temperature (T); therefore, with the 
adopted definition of the ionization pa- 
rameter, PocT/U. Under the assumption 
that the two components lie at the same 
distance from the central source, as is plau- 
sible given their identical kinematics, the 
two phases turn out to have the same 
pressure (see Table 3). To further illus- 
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Fig. 16. — Curve of thermal stability for 
the SED used in this analysis, as described 
in §4.1. Pressure equilibrium between the 
HIP and the LIP can be observed.^ 

trate this point, we show in Figure 16 
the thermal equilibrium curve (S curve) 
for the SED used in our analysis (see 
§4.1). To construct this curve we cal- 
culated the ionization parameter as de- 
fined by Krolik, McKee & Tarter (1981): 
S = [L / Airr"^) / c{nkTe) (where L is the lu- 
minosity of the source, r the distance be- 
tween the gas and the source, c the speed of 
light, n the total particle density of the gas, 
k the Boltzmann constant, and Tg the elec- 
tron temperature). This representation of 
the ionization parameter is simply the ra- 
tio between the radiation pressure and the 
gas pressure. In Figure 16 we also plot the 
S values we find for the HIP and the LIP 
(as listed in Table 3). The pressure equi- 
librium is evident. Unless imposed by real 
pressure equilibrium, this pressure equal- 
ity would have to be assumed to be a co- 



incidence. We notice that the HIP lies in 
the intermediate part of the curve, in a 
region where the slope is negative. Re- 
gions of the curve with negative slope are 
thermally unstable. This might be then 
inconsistent with the equilibrium scenario 
just presented. Nevertheless, the shape of 
the thermal equilibrium curve strongly de- 
pends on the SED, and particularly on the 
unobservable region between 100 A and 
912 A. Considering different shapes for the 
SED in this unknown gap can give rise to 
thermal curves with an intermediate sta- 
ble region in the region occupied by the 
HIP (see for instance Reynolds and Fabian 
1995). Therefore, this apparent instabil- 
ity might be only the result of our chosen 
SED, a likely possibility, since according 
to the error bars in Figure 16, our HIP is 
consistent with the semi-stable intermedi- 
ate region of our S curve (nearly vertical 
region of the equilibrium curve with a S 
value close to 1 and T 10^ K). Alter- 
natively, this instability could be refiect- 
ing the fact that, besides the illumination 
from the central source, additional cool- 
ing/heating processes (for instance shock- 
heating) might be also working on the gas. 
In this latter case, the electron tempera- 
ture would be constant due to the addi- 
tional source of heating, and the solution 
found would be stable. 

A two phase medium is consistent with 
the thickness ratio between the phases: 
the radial thickness occupied by the ab- 
sorbing gas is roughly equal to the ra- 
tio between the H column density and 
the H density (D=Ni://n(H)), and, from 
the definition of the ionization parameter 
given in §3, DcxNj:/r^U/Q(H). Therefore 
the thickness ratio between the phases 



33 



goes as Rtefcoc Nj^(HIP)r(HIP)2U(HIP)/ 
Niy(LIP)r(LIP)2U(LIP). If we assume again 
that both phases he at the same distance 
from the central source, then the thickness 
of the HIP is 140 times larger than that 
of the LIP. This result is consistent with 
pressure confinement of the LIP by the 
HIP. 

A two phase warm absorber with no in- 
termediate value of U was also detected in 
IRAS 13349+2438 by Sako et al. (2001). 
In that case, the UTA was also produced 
by Fe vil-xil. Are the two phases in IRAS 
13349+2438 also in pressure balance? An 
intriguing possibility is that a two phase 
medium forms because of the thermal in- 
stabilities. In photoionization equilibrium 
the gas will always be driven to the stable 
regions of the S curve which are set by the 
SED and atomic physics. However, not all 
the regions of the stable branches can be 
occupied by the absorbing gas, as the two 
phases will tend to reach pressure balance 
with each other. 

This appealingly simple picture not only 
suggests the idea of a two phase medium, 
but also points to a simple view of active 
nuclei (e.g. the structure for quasars sug- 
gested in Elvis 2000). Recent results in- 
dicate that the physical properties of the 
optical 'High Ionization Nuclear Emission 
Line Region' (HINER) are consistent with 
those of the X-ray warm absorber (Erkens, 
Appenzeller, & Wagner 1997; Nicastro et 
al. 1999; Nctzcr et al. 2002). Therefore, 
it is plausible that the high ionization ab- 
sorption and emission arise from the same 
gas. If the gas lies in our line of sight, the 
absorption is observed, if the gas lies out of 
our line of sight, only the coronal emission 
from the HINER can be detected. Por- 



quet, et al. (1999) found that this picture 
could only be supported if the electron den- 
sities are high (ug ~ 10^° cm~^) to avoid 
overproducing the coronal lines, and con- 
cluded that the distance of the ionized ab- 
sorber from the incident radiation source is 
of the order of that of the Broad Line Re- 
gion (BLR). All this evidence also points to 
models where the hotter gas confines the 
broad emission lines gas (BELs gas, see, 
e.g. Turner et al. 1993; Elvis 2000). 

Our two phase solution for the warm ab- 
sorber may also be pointing to this gen- 
eral scenario. The value of the electron 
density derived from the O VII triplet for 
NGC 3783 (see §4.2) is consistent with 
the predictions by Porquet et al. (1999, 
see above), and therefore this is consis- 
tent with a connection between the ab- 
sorbing and emitting gas. This density, as 
well as the temperature for the HIP are 
close to those needed to explain strong O 
VII emission lines in other objects (George, 
Turner, & Nctzcr 1995; Nicastro, Fiore, & 
Matt 1999). However, what is the loca- 
tion of the absorber? We do not have any 
means of measuring the electron density of 
the absorber or its distance to the central 
source at this time. Nevertheless, Nicas- 
tro, Fiore, & Matt (1999) and Netzer et al. 
(2002) have shown through spectral vari- 
ability studies in NGC 4051 and NGG 3516 
respectively, that the location of the ion- 
ized absorbing gas is similar or closer to 
the central source than that of the broad 
line region. 

A simple identification between the ab- 
sorbing gas and an emission component 
is difficult to obtain, but from the latter 
discussion a natural candidate could be 
the Broad Emission Line region. This has 
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been already suggested by Kuraszkiewicz 
& Green (2002), who found a correlation 
between the N v/C IV ratio in broad emis- 
sion lines and that in narrow absorption 
lines, and is further supported by the simi- 
larity of the temperature of the LIP and 
the temperature predicted for the BEL 
(e.g. Kaastra et al. 1995, for NGC 5548). 
The BEL gas in AGN is clearly radially 
stratified, as shown by reverberation map- 
ping (Peterson et al. 1999), with a range 
of ionization parameters. According to the 
above models, the BELs are produced by 
the contribution of all the clouds surround- 
ing the AGN, that can have different ion- 
ization conditions (simply because of dif- 
ferent distances to the central source, or 
different electron densities), giving rise to 
the range of ionization parameters. How- 
ever, the absorption is produced only by 
those few clouds intercepting our line of 
sight. Therefore, a discrete value of the 
ionization parameter for the LIP can in 
principle be reconciled with the LIP-BEL 
gas connection. Koratkar & MacAlpine 
(1992) found that a log(U) range between 
-1.3 and -0.7 could characterize the UV 
BEL in NGC 3783. This range could be 
consistent with the value obtained for the 
LIP. However, different SEDs were used in 
the studies, and though they have similar- 
ities (for instance, we included a change 
of slope in the UV, and they included a 
big blue bump) it is not possible to make 
a clear connection (see §5.4). Neverthe- 
less, the similarity of U (also found in other 
AGN, Elvis 2000), as well as the similarity 
of the temperatures are physically sugges- 
tive, and bear further investigation. 

This would point to models where the 
hotter emitting/ absorbing gas confines the 



BELR/ LIP (see, e.g. Turner et al. 1993, 
Elvis 2000). Pressure confinement is cur- 
rently not favored (Peterson 1997) among 
other reasons because of the rapid destruc- 
tion of the BEL clouds by drag forces, 
when they are moving through a stationary 
medium. However, pressure confinement is 
possible if both phases are comoving in a 
wind with a single outflow velocity, as is 
the case of our two phases on NGC 3783. 
Other problems, such as the need for a low 
Compton depth are solved if a funnel shape 
geometry is present (Elvis 2000). This ge- 
ometry is not in contradiction with our re- 
sults. 

As a final remark, we note that the pres- 
ence of the absorption component detected 
in O VII with an outflow velocity of 1365 
km s^^ is not straightforwardly explained 
in light of the above models. This com- 
ponent is well constrained, as discussed in 
section 4.3.2, although contributing to only 
4 features out of the 265 with EW larger 
than 1 mA. However, the thickness of the 
LIP is ~9 times that of this component, 
and that of the HIP is ~ 1200 times larger. 
Therefore, whatever the origin of this ab- 
sorption, it is not in contradiction with 
our scenario in which two phases domi- 
nate. Nevertheless, what is the true nature 
of this puzzling component? The velocity 
matches that of the component discovered 
in the UV, and in several other AGN mul- 
tiple velocity UV absorption systems have 
also been detected. Are these intervening 
systems unrelated to the AGN? Do they 
originate at the edge of an obscuring torus? 
The present data are inadequate to say. 

Our results, though strongly pointing to 
the general scenario presented above, are 
based only on the analysis of one object. 
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Was this result only a coincidence? Or 
can all ionized absorbers be represented 
with a simple two phase medium? IRAS 
13349+2438 is another example of an ion- 
ized absorber with two different absorp- 
tion components (Sako et al. 2001). Are 
these two components always in pressure 
balance? Further study is still necessary to 
unravel the answer to these questions, but 
the simple picture presented here is cer- 
tainly intriguing. 

6.2. The Nature of the UV X-ray 
Absorption 

A warm absorber was first discovered by 
Halpern in 1984, but the use of warm ab- 
sorbers to understand the physical condi- 
tions in matter close to the central black 
hole came with the realization that the 
warm absorbers seen in X-rays also imprint 
signatures in the UV as associated absorp- 
tion lines. In a series of papers Mathur and 
collaborators (1994 through 2001) showed 
that the high ionization absorption lines 
seen in the UV, such as Ovi, CiV, Lja 
could come from the same single phase gas 
which produced O Vll and O Vlll edges 
in X-rays. The combination of saturation 
free X-ray column densities with outflow 
velocities from UV spectra led to the dis- 
covery of a component of nuclear material 
not previously recognized: highly ionized, 
outflowing, low density material situated 
close to or outside the broad emission line 
region. The mass outflow rates implied are 
large, a significant fraction of the mass ac- 
cretion rate needed to power the AGNs, 
and so must have important dynamical in- 
fluences within the active nuclei. The in- 
sight obtained through these studies led to 
one of the most comprehensive models of 



the nuclear regions of AGNs to date (Elvis 
2000). 

As such, the X-ray/UV absorbers have 
proved to be powerful tools to understand 
the physics of AGN outflows. However, 
the connection between UV and X-ray ab- 
sorbers has remained a matter of contro- 
versy for years, for several reasons: 

(1) Intrinsic observational limitations 
make component identiflcations difficult: 
even now, UV spectra have 20 times bet- 
ter spectral resolution, so a simple, one to 
one, matching of features cannot be per- 
formed. In addition. X-ray instrumenta- 
tion is sensitive to only a limited range 
of column densities (see Fig. 11), making 
weak components found in the UV unde- 
tectable, even if present. For an ioniza- 
tion parameter U~ 1, FUSE can detect 
100 times smaller Nh systems (Kriss et al. 
2000) than even high S/N Chandra grat- 
ing spectra (Fig. 11). The situation gets 
more complicated because only a limited 
range of high ionization states can produce 
features in the UV spectra (Fig. 11), since 
the low ionization ion populations where 
the UV transitions are produced become 
tiny. We now see that a 2-phase medium 
is required. The high ionization phase in 
our model, as well as the high velocity sys- 
tem in NGC 4051 (CoUinge et al. 2001) 
represent good examples of this. 

(2) As already stressed by Mathur et al. 
(1994, see also Kaspi et al. 2001), the ac- 
tual ionization parameters are strongly de- 
pendent on the far UV shape of the SED. 
Although this does not affect the absorp- 
tion produced in the X-ray region (Steen- 
brugge et al 2003), it has an important ef- 
fect on the UV absorption and so on link- 
ing the UV and X-ray absorbers. This is 
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an important limitation since this region 
cannot be observed due to Galactic absorp- 
tion. 

(3) The strong variability of the sources 
makes the determination of the relative 
UV and X-ray continuum fluxes unreliable, 
unless truly simultaneous observations are 
considered (Crenshaw et al. 2003). 

(4) As discussed before, several studies 
on the UV have neglected the efi^ects of 
saturation, unc/erBStimating the ionic col- 
umn densities in the UV band. As demon- 
strated by Arav ct al. (2003), when an 
accurate treatment of saturation is present 
the column densities found in the UV band 
become consistent with the X-ray absorb- 
ing material. We flnd here that X-ray col- 
umn densities have instead overestimated 
the Oxygen edge depths and column den- 
sities due to line blends. 

(5) Finally, the presence of a hotter 
phase (not detectable in the UV), as well 
as the presence of deep bound-free absorp- 
tion edges in the X-rays, introduces ad- 
ditional complexity. While the UV ab- 
sorbers are almost transparent one to the 
other, photoelectric absorption in the X- 
ray attenuates signiflcantly the continuum 
each absorber receives. Hence, the order 
in which the absorbers are illuminated by 
the source may affect the final results in 
a given model. The ionization parameter 
is very sensitive to this effect, making the 
comparison of UV and X-rays values unre- 
liable. Therefore, to really infer something 
about the X-ray UV absorption link, simul- 
taneous self-consistent modeUing of each 
absorber should be applied to simultane- 
ously observed data. This is beyond the 
scope of this paper (but we notice that si- 
multaneous observations with FUSE and 



HST-STIS do exist for the Chandra data 
analyzed here). 

Some of the difficulties mentioned above 
can be overcome with careful analysis and 
modelling. But the complexity of UV and 
X-ray spectra have led to many papers 
with highly complicated models involving 
multiple physical components to describe 
AGN absorbers. The general consensus 
from these studies, using very long ob- 
servations on the most powerful satellite 
borne telescopes, seems to be that the ab- 
sorption system(s) is/are extremely com- 
plex. This renders warm absorber studies 
completely ineffective as probes of circum- 
nuclear environment of AGNs. 

In this paper we have shown that the sit- 
uation is much more promising: the warm 
absorber in NGC 3783 can be simply de- 
scribed by two gas zones, but no more, 
with identical kinematics within the power 
of our measurements (§4.3.1), one having 
a low and one having a high ionization pa- 
rameter; these two zones are in pressure 
equilibrium with each other and lie on the 
(T,S) equilibrium curve for the SED of 
NGC 3783, and so are two phases of the 
same medium; the LIP phase also produces 
the UV absorption lines (§4.4). 

Thus, the high resolution Chandra and 
XMM-Newton spectra can resolve years 
of controversy surrounding the nature 
of X-ray/UV absorbers, if analyzed and 
interpreted carefully. It appears clear 
that (though multiple components may be 
present) the X-ray and UV absorbers are 
different manifestations of the same out- 
flow (e.g. Steenbrugge et al. 2003), con- 
firming the original Mathur et al. models, 
but refining them to a two phase medium. 
This is extremely important, since the long 
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baseline of ionization states from UV to 
X-rays, and superior velocity information 
from higher resolution UV spectra, can be 
used to model the absorbing outflow accu- 
rately. We will exploit this tremendous po- 
tential in a forthcoming paper presenting 
detailed model of the X-ray/UV absorber 
in NGC 3783, and will extend the current 
analysis to more AGNs as high quality X- 
ray spectra become available. 

7. Summary 

We have developed PHASE, a new code 
designed to model the X-ray and UV ab- 
sorption of ionized gas in the environment 
of AGN. We have used our model to study 
the Chandra 900 ksec spectrum of NGC 
3783. This X-ray spectrum is the best one 
available so far from a Seyfert galaxy. We 
used a global fit approach, which is necessi- 
tated by the blending of 75% of the absorp- 
tion features, compromising line strengths 
measured empirically. The heavy blending 
also results in a continuum level that can- 
not be estimated locally, except in narrow 
bands. Our main results are as follows: 

(1) The intrinsic continuum of the 
source is well reproduced by a power law 
(r = 1.53) and a thermal component 
(kT=0.1 keV). We attenuated this con- 
tinuum by an equivalent hydrogen column 
density of 1.013 x 10^^ cm~^ to account for 
the Galactic absorption. 

(2) Our absorption model leads to a sim- 
ple picture consisting of a two phase wind. 
The model can reproduce more than 100 
features with only 6 parameters. The high 
quality of the data, plus the inclusion of 
new atomic data for several lines made the 
modelling results more reliable than previ- 
ous studies. The two components of the 



model are in pressure equilibrium, and are 
consistent with a single outfiow 750 km 
s~^), a single turbulent velocity (300 km 
s~^), and solar abundances. The ioniza- 
tion parameter of the high ionization phase 
(U=5.754) is 35 times larger than that of 
the low ionization phase (U=0.166). The 
equivalent H column densities were esti- 
mated as N// = 1.6 X 10^^ cm"^ for the 
hotter phase and Nj/ = 4.1 x 10^^ cm~^ 
for the cooler one. The difference found 
for the ionization parameters of the phases 
is strongly dependent on the continuum 
shape. With our chosen SED, a factor of 35 
yields a low ionization phase with a tem- 
perature 37 times lower than the high ion- 
ization one, and completely different ion- 
ization degrees. The main features of the 
cooler phase are the Fe M-shell UTA and 
the O VII lines, while the features for the 
hotter phase are the O Vlll and the Fe L- 
shell lines. The O VII, previously identified 
with the O VIII and a hot phase, is consis- 
tent with a cooler phase and the O VI. 

(3) Our model strongly disfavors a con- 
tinuous range of ionization parameters to 
describe the medium because it is not con- 
sistent with the observed shape of the 
UTA. 

(4) Thanks to the high resolution spec- 
tra, we were able to determine that in the 
case of NGC 3783, the K edge of O vii 
is masked by the presence of the UTA and 
the K edge of O VIII is masked by the pres- 
ence of the Fe L-shcll lines. Through our 
model, we have shown that this produces 
an overestimation in the measurement of 
the edge depths (by a factor > 4), if the 
contribution from the absorption lines is 
not taken into account. An important con- 
clusion derived from this is that only an 
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upper limit to the column densities can be 

obtained from CCD low resolution spectra 
(e.g. ASCA). In objects where a deep UTA 
is present, this effect could be particularly 
important. 

(5) Our model predicts a Ca xvi hne at 
21.45 A. This line was formerly identified 
as a local O vil(A21.602) feature arising 
from an intergalactic cloud at zero rcdshift 
or Galactic absorption. With the inclusion 
of the Ca hne, if a local feature is present, 
it cannot contribute more than 9.3 mA to 
the total EW. 

(6) Our model does not attempt to pre- 
dict the emission lines in a self-consistent 
way. Therefore, we fitted them with gaus- 
sian profiles. We restricted the outfiow 
velocities and FWHMs of the lines aris- 
ing from transitions of the same ion. The 
emission and absorption lines are blended, 
making an independent estimation of them 
unreliable. Our emission lines are con- 
strained by the absorption lines, making 
our predictions useful as rough estimates of 
the actual emission line properties of NGC 
3783. The emission line EWs increase by 
a factor ~2 when the absorption/emission 
filling effect is included. 

(7) We find an excellent agreement be- 
tween X-ray data and a model with kine- 
matic components identical to the UV ab- 
sorption hnes. We also find that the UV 
value of ionization parameter underpre- 
dicts ionic column densities by as much 

factor of ten, illustrating the danger 
of modelling the absorbers using UV data 
alone. A detailed model of the UV and X- 
ray data together, however, is beyond the 
scope of this paper. 

The analysis presented in this paper of- 
fers new hope that warm absorbers in AGN 



are relatively simple systems with physical 
properties that are easily understood. If 
this turn out to be generally true, then 
they will provide, particularly through 
their variability, completely characterized 
physical conditions for AGN winds. This 
in turn is likely to lead back to a deeper un- 
derstanding of the basic physics powering 
AGN and quasars. 
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